HAUSDORFF MEASURE OF VORTICITY NODAL SETS FOR 
THE 3D HYPERVISCOUS NAVIER STOKES EQUATIONS WITH 

GENERAL FORCES 

ABDELHAFID YOUNSI 



Abstract. In this paper, we modified the three dimensional Navier-Stokes 
equations by adding a i-Laplacian. We provide upper bounds on the two- 
dimensional HausdorfT measure Hf of N® = {x G £7 C K 3 / uj(x,t) = 0} the 
level sets of the vorticity lj of solutions. We express them in terms of the 
Kolmogorov length-scale and the Landau-Lifschitz estimates of the number 
of degrees of freedom in turbulent flow. We also, under certain hypothesis 
recover the two-dimensional HausdorfT measure estimates for the usual 3D 
Navier-Stokes equations with potential force. Moreover, we show that the 
estimates depend on I, this result suggests that the modified Navier Stokes 
system is successful model of turbulence and the size of the nodal set Hf(N®) 
leads the way for developing the turbulence theory. 



1. Introduction 

In this paper, we provide upper bounds on the two-dimensional HausdorfT mea- 
sure Ttf of N® the level sets associated with the vorticity of modified three di- 
mensional Navier-Stokes equations. We modified the 3D Navier-Stokes system by 
adding a higher-order viscosity term to the conventional system 

^ + e (- A)' u - vAu + (it.V) u + Vp = f (x) , in fi x (0, oo) 

divw = 0, in x (0, oo) , h .1) 

p(x + Lei, t) = p(x, t), u(x + Lei, t) — u(x, t) i — 1, d t S (0, oo) 
u (x, 0) = uq (x) , in n, 

on ft = (0,L) d with periodic boundary conditions and (e\,...,ed) is the natural 
basis of M. d . Here e > is the artificial dissipation parameter and v > is the 
kinematic viscosity of the fluid, I > 1. The functions u is the velocity vector field, 
p is the pressure, and / is a given force field. For e — 0, the model is reduced to 
the Navier-Stokes system. 

In the work [31] , the strong convergence of the solution of this problem to the 
solution of the conventional system as the regularization parameter goes to zero, 
was established for each dimension d < 4. 

Mathematical model for such fluid motion has been used extensively in turbu- 
lence simulations (see e.g. [5]) also see Borue and Orsag [31 2]. For further discussion 
of theoretical results concerning (11.11 ). see [II l2l [23l [31] . 

For the 3D Navier-Stokes system weak solutions of problem are known to exist 
by a basic result by J. Leray from 1934 [21j . only the uniqueness of weak solutions 
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remains as an open problem. Then the known theory of global attractors of infinite 
dimensional dynamical systems is not applicable to the 3D Navier-Stokes system. 

In particular, in case one accepts the point of view that the dimension of a global 
attractor for the Navier-Stokes equations is associated with the number of degrees 
of freedom in turbulent flows, then the two-dimensional Hausdorff measure Hf(N^) 
is an important way to the understanding of turbulence theory [31] . 

We are interested in the three dimensional case. Let P m be the projection onto 
the first m eigenspaces of the Stokes operator A = — A and let N® = {x £ ft C K 3 / 
u>(x,t) = 0} the nodal sets of the vorticity u for solutions of the equation (|1.1| ). 
We provide an upper bound on the size of the nodal sets Uf{N^) and we show that, 
the bounds necessarily depend on m and I this dependence is a fractional power 
of I. Thus answering a question raised by J. Avrin pQ. We also obtain here scale- 
invariant estimates on the two-dimensional Hausdorff measure Hf(N%) in terms 
of the Landau-Lifschitz theory of the number of degrees of freedom in turbulent 
flow. Since expressing the above estimates in terms of the (dimensionless) Grashoff 
number G. In order to obtain an upper bound on the Hausdorff measure of level 
sets associated, we use the method from [19] (see also [6], [7]). 

The main purpose of the present article is to study the dependence of the two- 
dimensional Hausdorff measure Hf{N®) on the parameter I. Using a family of 
Kolmogorov flows as base flows we can deduce also upper bounds on the Hausdorff 
measure Hf(N®). We also find here that the upper bounds on the two-dimensional 
Hausdorff measure of N° converges to the corresponding upper bounds on W\{N®) 
the two-dimensional Hausdorff measure of the nodal sets of the usual 3D Navier- 
Stokes as I — 1. Under certain hypothesis we recover the two-dimensional Hausdorff 
measure H\(N®) estimates for the usual 3D Navier-Stokes equations with potential 
force. We extend the method from [17] to a 3D Navier-Stokes with general forcing 
modified by Z-Laplacian. These estimates are obtained without using the Dirichlet 
quotients [IT) . 

We note, however, that for the incompressible 3D Navier-Stokes equations with 
general force, it seems not so easy to get some better estimates on the Hausdorff 
measure of the level sets associated with the vorticity as in the case of potential 
force studied in [TTJ [18] for periodic solutions of the 2D. Related results for the 
3D Navier-Stokes equations (with general forcing) can be found in !5]. The upper- 
bounds on the Hausdorff measures of the level sets associated with solutions of 
some other partial differential equations were obtained in [TJ,[TT], [T^], [IB], [19] , 
[20), and [22]. 

The paper is organized as follows. In Section 2, we present the relevant math- 
ematical framework for the paper. In Section 3, we provide upper bounds for the 
two-dimensional Hausdorff measure Hf of the level sets associated with the vorticity 
of the Navier-Stokes system with hyperdissipation. 



2. Notations and preliminaries 

In this section we introduce notations and the definitions of standard functional 
spaces that will be used throughout the paper. We denote by H m (O), the Sobolev 
space of Lper periodic functions. These spaces are endowed with the inner product 

(u,v)= (£ )/3 u,D' 3 w) L 2 (0) and the norm ||w|| m = £ (||-D^|| i2(n j5. 

\P\<m \P\<m 
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Each u G L per can be identified with its Fourier expansion 

u (x) = Uk cxp(2ink.—) 

kez, 3 

where G C 3 satisfy Uf, = Then u is in L 2 if and only if 

h* 3 = m £ i«*i a < «>, im = x 3 , 

fcGZ 3 

then the Sobolev space u £ H m (17), m G K + can be characterized by 
ff m (17) = { M ,H fe = ^ fc 2m |u fc | 2 < oo.} 

H- ,n (17) denote the dual space of Zf m (17). 

We denote by i? m (17) the subspace of H m (17) with, zero average 

H m (17) = {u G H m (17) ; Ju(x)dx = 0}. 

n 

For m = 0, we have £T m (17) = L 2 (17). 

• We introduce the following solenoidal subspaces V s , s £ R + which are 
important to our analysis 

V (17) — {u £ L 2 (17) ,divu = Q,u.n | E< = -u.n \s i+3 ,i = 1,2,3}; 

Vi (fi) = {u £ H 1 (17) , divtt = 0, 7o u | Ei = 7o^ |s I+3 , i = 1, 2, 3}. 

V2 (H) = G H 2 (17) ,divu = 0,7 u |e<= 7o^ |s i+3 ,7i w |e«= ~7i w |s i+3I * = 1,2,3}, 

see [29l Chapter III, Section 2]. We refer the reader to R.Temam [30] for details 
on these spaces. Here the faces of 17 are numbered as 

Si = dQ. n {2^ = 0} and S l+3 = <9ft n {x, = L} , i = 1, 2, 3. 

Here 70, 71 are the trace operators and n is the unit outward normal on 917. 

• The space Vq is endowed with the inner product (it, i>)z, 3 (n) an< ^ norm 
INI = (w.«)i?(n)- 

• Vi is the Hilbert space with the norm HmHj = ||u||y ■ The norm induced by 
H 1 (17) and the norm ||Vu|| are equivalent in V\. 

• V2 is the Hilbert space with the norm = \\u\\ v ■ I n V2 the norm 
induced by H 2 (17) is equivalent to the norm ||Au||. 

V' a denote the dual space of V s . 

Let P be the orthogonal projection in L 2 cr (M 3 ) 3 with the range H. 
Let A = —PA the Stokes operator. It is easy to check that Au = — Au for every 
u £ D (A). We recall that the operator A is a closed positive self-adjoint unbounded 
operator, with D (A) = {u £ Vq, Au £ Vq}. We have in fact, 

D[A) = H 2 (17) n V = V 2 . 

The eigenvalues of A are {Aj}^~^, < Ai < A2 < ...and the corresponding or- 
thonormal set of eigenfunctions {wj}^^ is complete in Vq 

Awj = XjWj, Wj £ D(A),\/j. 
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The spectral theory of A allows us to define the powers A 1 of A for / > 1, A 1 is an 
unbounded self-adjoint operator in Vq with a domain D(A l ) dense in V2 C Vq. We 
set here 

A l u = {-A) 1 u for u (z D (A 1 ) = V 2l n V . 
The space D (A') is endowed with the scalar product and the norm 

( u > v )d(a>) = (A l u,A l v) , \\u\\ D(Al) = {(u,w) D(A1) }i (2.1) 

00 

In the case for I > 0, we have D [A 1 ) = {u € g,y^A^ (tt,^) 2 < 00}. For / e K 

j=i 

the scalar product and the norm in (|2.1I ) can wiriten alterntivly as 

00 00 
{u,v) D{Al) =^2xf(u,w j ){v,w j ),\\u\\ D{A , ) ={^Af (2.2) 
3=1 3=1 

and for u € D(A') we can write 

00 

A'u = /~^A j (it, Wj)wj. 
3=1 

Let us now define the trilinear form &(., ., .) associated with the inertia terms 

3 

b(u,v,w) = ^ 1 



n 



The continuity property of the trilinear form enables us to define (using Riesz 

representation Theorem) a bilinear continuous operator B (u,v); V2 x V2 — > will 
be defined by 

(B(u,v),w) = b(u,v,w), G V2. (2.3) 

Recall that for u satisfying V.it = we have 

6 (it, u, u) = and 6 (u, u, w) = —b (u, w, v) . (2-4) 

We recall some well known inequalities that we will be using in what follows. 
Young's inequality 

ab < -aP + -^-rb q ,a, b, a > 0, p > 1, q = (2.5) 
P qcrp ' P-1 

Poincare's inequality 

Ai ||u|| 2 < \\A^u\\ 2 for all u E V . (2.6) 

Denoting 

NlG(t) =11 ^ II and (u,v) G(t) = (e tA K,e tAi v). 

The set D(e aA ) is called the Gevrey class of operator of order a > [TU]. Our use 
of Gevrey classes shall be based on the following consideration. 

Denote with N® = {x e £1 : h(x) = 0} the zero (nodal) set of a function ft in a 
set f2, and let Ti 2 be the two-dimensional Hausdorff measure operating on subsets 
of R 3 (area in this case) . 
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3. Level Sets of the Vorticity Function 

Using the operators defined above, we can write the modified system ljl.ll ) in 
the evolution form 

du 

— + eA'u + vAu + B (u,u) — f (x) , in SI x (0, oo) ,^ ^ 

u a ( x ) = wo, in SI. 

The existence and uniqueness results for initial value problem (11.11 ) can be found 
in |23j . The following theorem collects the main result in this work 

Theorem 3.1. Fori > d is the space dimension, for e > fixed, f € L 2 (0, T; Vq) 
and uq £ V$ be given. There exists a unique weak solution of ) which satisfies 
u e L 2 (0, T; Vi) n L°° (0, T: V ) , VT > 0. 

The modern understanding of turbulence is that it is a collection of weakly cor- 
related vortical motions, which, despite their intermittent and chaotic distribution 
over a wide range of space and time scales, actually consist of local characteristic 
'eddy' patterns that persist as they move around under the influences of their own 
and other eddies' vorticity fields [15] . 

In fluid mechanics, the Reynolds number is important in analyzing any type of 
flow when there is substantial velocity gradient (i.e. shear.) It indicates the relative 
significance of the viscous effect compared to the inertia effect. The Reynolds 
number is proportional to inertial force divided by viscous force (see [9] ) 

Re=^- U 2 = L- 2 (\\u\\l) (3.2) 

where I the characteristic scale of the forcing and (.) is the long-time-average 

(<?(.)) = r lim sup(i fig(t)dt). (3.3) 

With Reynolds number calculator we can analyze what makes fluid flow regime 
laminar and what is needed to force the fluid to flow in turbulent regime. Experi- 
mental observations show that for 'fully developed' flow, laminar flow occurs when 
Re < R t e and turbulent flow occurs when Re > R t e. In the interval between R t e 
and R t e, laminar and turbulent flows are possible ('transition' flows) |9i and refer- 
ences therein. The nature of the vortex formed in the fluid flow depends strongly 
on the Reynolds number([5]; and references therein). These transition Reynolds 
numbers are also called critical Reynolds numbers, and were studied by Osborne 
Reynolds around 1895 [55]. The transition to turbulence and the constructon of 
vortex are delayed by increasing the critical Reynolds number. If we assume that 
the critical Reynolds number R c e for the onset of vortex shedding is, atteint for 

11-11 = ^L, (3.4) 
then the associet velocity u for each 

Re > R c e (3.5) 

satisfaies the inequality 



\\u\\>^L = n, (3.6) 



is a positive constant. 
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Another nondimensional quantity that we use often is the so-called Grashof 
number, which is proportional to the forcing term /. Hence, we define the Grashof 
numbers in the 3-dimensional case, as in Foias, Manley, Rosa and Temam [5] by 

Gr(/) = ^7574 11/11 (3 ' 7) 

The effects of variation in Grashof number on vortex have been shown in the work of 
Olson and Titi 23] , they keep the spatial structure of the forcing function fixed and 
vary the Grashof number by varying the amplitude of the forcing function. Namely, 
they vary the Grashof number by rescaling the forcing function by a multiplicative 
factor. This is equivalent to changing the viscosity or the size of the domain. 
As increases, or equivalently as the viscosity decreases, the turbulent flow becomes 
more energetic and one would expect the number of numerically determining modes 
to increase as well. There are many reasons to suppose that the existene and 
intensty of vortex in our work should increase as the grashof numbre increases 
[THIS, 25 . In [M] zero forcing implies that the attractor has been reduced to zero. 
Since all solutions decay eventually to zero in the unforced case. 

This intuition is supported by existing theoretical critucal numbre G c r(f) for 
the existence of level curves of representative vorticity fields. 

Note the flow for Gr(f) > G c r(f) has noticeably more large scale structure 
compared to the flow for Gr(f) < G c r(f). This is consistent with the energy 
spectra, where most of the energy is in the lowest modes, that is, in the large 
spatial scales and eddies when the Grashof number is large [24| ■ 

The effect of a body force on vorticity production and turbulence generation in 
a fluid flow is described by the Grashof numbre. 

In addition, we assume without loss of generality that ||/|| is bounded. Than, 
there exist a maximum Grashof numbre G max r(/) and a positive constant p such 
that the body force / satisfies the follwing inequality 

11/11 <^ 2 A? /4 G max r(/)=p. (3.8) 

Since ||/|| is srictement positive we get 

11/11 P ^A? /4 ZG max r(/) 
this gives a relation between ||u|| and ||/|| 

\H>P 11/11- (3-io) 

Moreover, according to the definition of the Gevrey norm and the relation (|3. 101 ) 
we get 

NI G >/3||/II G - 
The vorticity, uj = V x u satisfies the equation 

(4 + u.V + z/A + e (-A) l )u = lj.Vu + F (3.11) 
dt 

where F = V X /. 

Theorem 3.2. [19] Suppose that a nonzero function h G V\ satisfies 

We^hW^MWhW, 

Then 

H 2 (iVg) < C1L 2 (1 + log Af) e C2L/a . 
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Hereafter, d for i € N, stand for universal constants. The above statement is a 
special case of [HI Theorem 2.1]. It will be used in conjunction with the following 
statement: 

Lemma 3.3. \17\Let u € Vq, and let uj and be its vorticity. If 

|| Ah aA K \\< M || A?u || (3.12) 
for some M > 0, then, for every c G M 

||e^*(u,-c)||<M||a;-c||. (3.13) 

For the rest of the paper, let u(t) be an arbitrary solution of the the modified 
Navier Stokes system (jl.ll ) with u(0) = uq. 

Theorem 3.4. Let \\u\\ > f3 \\f\\ for any a < — — and j3 < — — . TTien £/iere exists 

4 z/ 

a universal constant C3 such that if \\ A^u |j< C3, iften 

|| i5e aal !i ||< 2 || A?u ||, i > 0. (3.14) 
Proof. For any a, t > 0, We take the inner product of (|3.1I ) with w, to obtain 

\j t \\ A ^ u \\ 2 G(t) = a \\ A ~ iu \\ 2 G{ t) + {Au,u) G{t) 

= <4 A *u\\ G (t) - eP^^llcw - v \\ Au \\ G (t) - b ( u > u, Au) G[t) + (/, Am) 



then using the Young's inequality (|2.5I ) we have 



(3.15) 



^ll^ u ll 2 G( t ) < - £ \\A^u\\ 2 G{t) + ~\\Au\\ 2 G{t) + ^\\Aiu\\ 2 G(f) - v\\Au\\ G{t) 
+Q\M\ G ( t ) + ^WfWUt)) + b(u,u,Au) G(t) 



From we get 

this give 
Id,,,! 



\\M\hit) > ll/llcw 



2^ii— "HG( t) < -4A^u\\ 2 G{t) + ^\\Au\\ 2 G{t) + ^\\Aiu\\ 2 G{t) - v\\Au\\ 2 Git) 
+(^\\M\ 2 G (t) + ^\\^u\\ 2 G{t) ) + b{u,u,Au) G (t) 
< -e\\A^u\\ 2 G{t) + zhL\\ A iu\\ 2 G(t) + ^-||A*u||» {t) 



P 

2X^^ U ^G(t) +b{u,u,Au) G{t) . 

1 



We get for /3 2 



\l t \\ Ahu \\l {t) < -4 A ^A 2 m - ^f\\ A H 2 m 

2a 2 1 

H \\A2u\\ 2 G{t) + b(u,u,Au) G{t) . 



- .).-(/ 

v 

We use the following inequality from [10] and jTJl Section 4] 

\\Au\\ 2 

b(u,u,Au) G(t) <C4A^u\\ 2 G{t) \\Au\\ G{t) (l + log - 7' )° (3.16) 
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to obtain 



+ C 4 ||^.||^||A M || G(t) (l + lo gAi|| ^ u|| ^, 

To establish (|3.14l ) we use the estimate [17] 

a A i(l+log^ 7 )5 <d A i 2 + — log— a,d>0,^>6>0. (3.17) 
tr or ba 

By applying the Poincare's inequality (12.61 ). we have that for /i = || ^4u|| G(t) an d 
d 



|G(t) 
V 



IjtW^Wh) + 4A^u\\% (t) < -^\\Aiu\\% {t) + ^\\Aiu\\% (t) 

+C 5 \\A^u\\% (t) log Ce]lA h llG<t) . 

. v\\ , „ n 4^/2 , 

Letting a < — — we have for p < that 

4 v 

i|||^ U || 2 G(t) <C7 5 ||A^|| 4 G(t) log ^SS±. (3.18) 

i 

If || ^4 2 tt 1| < -4- = C3, the term with a logarithm in p. 181 ) is negative at t = 0, 
and thus (13.181 ) implies that ||^43u|| G /^ is a decreasing function of t. □ 

Theorem 3.4. implies that, for any solution u(t) of (jl.ll ). the space analyticity 
radius of u{t) goes to infinity as t — > 00. 

Let f2 be a periodic box, for simplicity assume fi = (0,L) 3 , A has eigenvalues 
< Ai < A 2 < ... with corresponding eigenspaces Ex, E 2 , ... Let P m be the 
projection on the eigenspaces E 1 (BE 2 (B---(BE rn and let Q m — I—P m we have||u|| 2 = 
||fmu|| 2 + ||Qmu|| 2 and we also have from (|2.2I ) that 

||A z n|| < X l m \\u\\ for every I > and u £ D{A l ). (3.19) 

For any t > 0, let u(t) be the vorticity of u(t). We shall, for any fixed t > 0, 
estimate the quantity 

I (w (t)) = sup H 2 (A^). (3.20) 

ceR 

Recall that for a function /i : O ->• K, N% = {x £ fi : h(x) = 0}. We need the 
following fact 

Lemma 3.5. Let \\u\\ > /3 ||/|| for any (3 > 0. Then 

\\u(t) || > ||ti(0)|| exp(?7t) for every t > 0. (3.21) 

With r] = -{exi + ^f-). 

Proof. Taking the scalar product of both sides of (jl.ll ) by u{t) and using (|2.4I ). 
we have that 

~\\u\\ 2 + v\\A?u\\ 2 + e\\A?u\\ 2 = (f,u) for t > 0. (3.22) 
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Using p. 191 ) and the following inequality 

(/,«)> -\\\f\\ 2 -\\\u\\ 2 (3.23) 
Because the increasing sequence < Ai < A < A m we have 

\j t \\nf > (-^A||A^|| 2 - (e\i + ^^)hl| 2 (3-24) 

note that since 



II A*u\\<C 3 

we have that 

lj t h\\ 2 >-c 3 -(exi + ^l)\\v\\ < : ^ : " 



2/3 

1 d 



if we set n = — (eX? n H 2 ) then we have from (|3.25l ) that 



Integrating the above inequality from to t, we get 

|2 ^ ~ C 3 ft „™/W\^ , IL./AMI2. 



H> (1 - expfat)) + ||«(0)|| 2 exp{r]t) (3.26) 

V 



or, since 

-c 3 



1 - exp^t)) > 0. 



Thus, we have the inequality (|3.21l ). □ 

i_ 

Proposition 3.6. Let \\u\\ > /3\\f\\ andu Q ^ 0, and suppose that \\uq\\ < C^vX^ .Then 
we have that 

-Lo 5 - + ( £ A^ + ^ • 



I (w (t)) < CiL(l + -Lo,?^ + (eA^ + _^)i) e — /or i > 0, (3.27) 



vX\ , „ 4^ 

for any a < ana a < . 

4 z/ 

Proof. By Theorem 3.4 and Lemma 3.5, we get for t > the following 

|| A3e atA *w || < 2 || Aiu || 
< 2A5 Huoll 

and use the inequality (|3.21l ) to get 



(3.28) 



(3.29) 



|| Ah atA ^u || < 2Xi \\ u (t)\\ exp(eA* + l -±f)t 

< 2(£)* || A*u(t) || cxp(eAi + ^f-)t. 

The rest follows by combining (|3.29l ) with Lemma 3.3 and Theorem 3.2. □ 



The foundational result for our two-dimensional Hausdorff measure estimates of 
= {x G £1 C M 3 / w(a;, f ) = 0} the level sets of the vorticity U) of solutions is 
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Theorem 3.7. Let \\u\\ > /3\\f\\ and uq ^ 0, and suppose that \\uq\\ < C31/X' . 
Then 

I (w (t)) < C 7 X? n for t>to, (3.30) 

2f , r f A? 4^2 

tuif/j io = 1 for cmw a < and 8 < . 

uxl 4 v 

Proof. With t > 2 ^ the inequality ([5~27l ) implies 

i/A? 



I (w (*)) < Cie(2 + -io 5 ^ + (eA^ + ~^-) — ) for t > 0. (3.31) 



-Lo. - + ( £ A^ + ^) — 

Since A m > Ai ()3.31l ) follows directly from the above inequality. □ 

The estimate of the Hausdorff measure Hf grows in m due to the term ^ but 
at a rate lower than | . 

Proposition 3.8. Let \\u\\ > f3\\f\\ and uo 7^ 0, and suppose that \\uq\\ < C^vX^ . 
Then 

sup I (w (*)) < Cgrna /or f > t . (3.32) 

£— >oo 

/or an?/ a < and p < . 

4 1/ 

Proof. Note that in the 3D case we have Xj > C$L~ 2 ji for some positive universal 
constant (see, for example j^HJ Lemma VI 2.1]). Therefore, Since X m ^ Aims the 
growth in m of the Hausdorff measure (|3.32l ) is less than mh . □ 

If we impose the condition A m < (j-) 2 or ^ < (j 2 -) 2 where io = A^ represents 

characteristic macroscopic length, and l e is the Kolmogorov length scale, i.e. l e = — 
where e is Kolmogorov's mean rate of dissipation of energy in turbulent flow (see 
e.g. [TJ IS EH US] , and the references contained therein) is defined as 

T 

e = Af f lim sup / || A* || 2 ds. 

T^co J 


Substituting this in (13.301 ) gives 

bu P Z(w(*))<Cio(^)*. (3.33) 

it— >oo 'e 

Since the (dimensionless) Grashoff number G = sup 'So|/ll2 in 3D ( gee e g [Jj[7jj[2g]) 

is an upper bound for { l f-) 2 . Hence, we obtain for the Hausdorff measure of the 
equation (jl.ll ) the following estimate in terms of the Grashoff number G. 

1 

Proposition 3.9. Let \\u\\ > j3\\f\\ and uq ^ 0, and suppose that \\uq\\ < C^vXl . 
1 

Then for any a < — — and 8 < — = we have 

sup I (u (t)) < C U G* for t > t . (3.34) 

t— >QO 
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This result holds independently of to, with C\\ independent of to. The estimate 

grows in I at a rate lower than |. If we impose the condition 1 = 1, the estimates 

i 

become sup^^ I {ui (t)) < CnG v . This result recover the usual 3D Navier-Stokes 
equations estimates, for the two-dimensional Hausdorff measure Wi(N®) estimates 
of the level sets associated with the vorticity. Here again our results indicate that 
under certain conditions the upper bounds for Hi{N^) converge to the associated 
upper bounds of the two-dimensional Hausdorff measure Hf(N°) estimates for the 
usual 3D Navier-Stokes equations with potential force. 

4. Conclusion 

Proving global regularity for the 3D Navier-Stokes equations is one of the most 
challenging outstanding problems in nonlinear analysis. The main difficulty in es- 
tablishing this result lies in controlling certain norms of vorticity. More specifically, 
the vorticity stretching term in the 3D vorticity equation forms the main obstacle 
to achieving this control, C. Foias [8] and estimates on the number of degrees of 
freedom for the Navier-Stokes equations and its closure models are a measure of 
the complexity of the system J. Avrin [T]. This paper proposed another interesting 
way to estimate this complexity through bounding the size of the nodal set for the 
vorticity and expressing this estimate in terms of G. 

We provide upper bounds for the two-dimensional Hausdorff measure TLf of the 
level sets associated with the vorticity of modified three dimensional Navier-Stokes 
equations this bounds depend on to and /, this dependence is a fractional power of 
I. Thus answering a question raised by J. Avrin pQ. 

Another interesting way to study decaying turbulence in the three-dimensional 
incompressible Navier-Stokes equations is to prvide a numerical investigation of our 
theoretical results on the size of the nodal set for the vorticity in the dependence 
of turbulence structure and vortex dynamics, as was done in |24j for the number 
of numerically determining modes in the 2D Navier-Stokes equations. It would be 
interesting to see how the turbulence structure depend on I. 
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